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Abstract. The magnetic- field dependence of many-body states in quantum dots 
can be tailored by controlling the mixing of various angular momenta. In lateral 
quantum dots — defined electrostatically in a two-dimensional electron gas this mixing 
can be accomplished by introducing anisotropics in the confinement potential, thereby 
explicitly breaking rotational symmetry. Mixing can be severe enough to violate Hund's 
rules, even at zero magnetic field. We illustrate the principle through calculations 
of states and spectra of four-electron droplets (p-shell) with long-range Coulomb 
repulsions and confined in anisotropic potentials. Our results show that the Hilbert 
space in these nanostructures can be engineered to the particular application domain. 



PACS numbers: 73.21.La, 73.23.Hk, 73.63.Kv, 85.35.Gv 
1. Introduction 

In order for quantum nanoelectronics to become a demonstrably viable technology, the 
particular nature of the discrete states of a finite system — a quantum dot, for example — 
must be both known and controlled. This is also a natural prerequisite for controlling 
the coherent time evolution of a nanoscopic quantum system. 

In the present work, we study the specifics of particular quantum dots jT] defined 
electrostatically within a two-dimensional electron system (lateral quantum dots). By 
adjusting the shape of the confinement potential, individual many-body states within the 
quantum dots, and within the leads, may be tailored to exhibit a desirable magnetic-field 
dependence. (See, for example, Ref. 0.) Among similar lateral structures, the present 
dots dots have the distinguishing feature that the number of electrons can be controlled 
down to a single electron pQ, similar to vertical devices 011]. They differ from vertical 
devices in that the confinement potential can be generally controlled to a greater degree; 
in the dots we study here, confinement is produced by four different gates. Because of 
this, singlet-triplet transitions [3] and negative differential resistance [H], for example, 
can be tuned at a fixed magnetic field by adjusting the gate voltage alone. Further, 
they differ from vertical devices in that the leads of lateral devices are the edges of a 

t URL: http://soliton.phys.dal.ca 
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high-mobility two-dimensional electron gas (2DEG). In such devices, spin-splitting of 
the edge states, due primarily to exchange effects, results in a spatial separation of the 
spin channels. Consequently, the tunnelling rates into and out of the quantum dot are 
measurably different for each spin species. This additional diagnostic tool has been 
recently used |Zj to show how interactions force the collapse of the spin-singlet state at 
filling-factor v — 2. 

The present lateral dots have a relatively small confinement energy (ujq = 1 meV) 
compared to, say, vertical devices. Because of this, the system is more strongly 

1 /2 

interacting (Coulomb interactions ~ u Q , kinetic energy ~ u ). This enhanced 
interaction strength largely drives the various phases we describe below. The dots are 
also anisotropic [HI EH EH and, consequently, the many-body states of the system contain 
mixed angular momenta. By studying the magnetic-field evolution of these states, their 
angular-momentum composition can be deduced. 

The dots are well-isolated from the leads and we work in the Coulomb blockade 
regime; transport proceeds via single-particle sequential tunnelling into and out of the 
quantum dot. We thus assume an isolated dot and look at its intrinsic properties. 

2. The isotropic regime 

Insight can be gained by looking first at the simple case of parabolic confinement. 
We first review the relevant spectra for dots containing four electrons. In these two- 
dimensional dots, the s-shell has one orbital with angular momentum L z = and the 
p-shell has two degenerate orbitals with angular-momentum L z — ±1. Application of a 
magnetic field perpendicular to the 2DEG breaks the orbital degeneracy. 
The noninteracting Hamiltonian is 



where r = (x, y) is the two-dimensional position operator. The single particle states are 
the well known fH\ IT2*] Fock-Darwin states, or alternatively , and more useful in the 
present context, a set of two decoupled harmonic oscillators \mna): 



In this notation, the ^-component of angular momentum of the state \mna) is L z = 
h(n — m), n is the Landau-level index, and 



with uj c = eB/(m*c) the cyclotron frequency and u the parabolic confinement 
frequency. We also take the magnetic field B in (JTJ), and throughout this paper, to 
be directed perpendicular to the two-dimensional plane of the dot. The final term in 
(J2J) is the Zeeman energy. 

When more than one electron is trapped in the dot, long-range Coulomb repulsion 
plays a significant role, particularly for few-electron droplets, where screening is 
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Figure 1. Ground-state transitions for 2D parabolic dots containing four electrons 
with confinement ujq = 1 meV. The figure shows those low-lying states which will be 
significantly mixed by the anisotropy. (See text.) The states shown all have S z = 
and material parameters are for GaAs. 

minimal and the long-range nature of the interaction is manifest. We take the full 
Coulomb interaction Vq = e 2 /(e|r 1 — r 2 |) and solve the problem numerically by 
exact diagonalisation, exploiting the spin and rotational invariance of the model. (See 
Refs. for details.) 

Figure [T] shows the successive ground states for four particles jS], for GaAs 
(material parameters: dielectric constant e = 12.4; effective mass m* = 0.067; g- 
factor g = —0.44). At zero field, the ground state configuration is a spin triplet 
with zero angular momentum. We label this state "Hund's-rule state" because the 
dominant component, approximately 75%, of this correlated state is the Hund's-rule 
state 1 00 |,00 f, 10 X, 01 |) §. The first crossing occurs already below 0.1 T, and 
the new ground state is a spin singlet state with angular momentum L z = —2. The 
dominant component of this state is the usual v = 2 state, 1 00 j, 00 t, 10 |, 10 f )• Hence, 
this transition involves both a spin and an orbital transition. The second ground-state 
transition of relevance here, occurring at approximately 0.45 T, is to an "edge triplet" 
state. This is again a spin and orbital transition; the new ground state is a spin triplet 
and has angular momentum L z = —3. The dominant contribution to this state is the 
single Slater determinant 1 00 J.,00 t, 10 j.,20 j). We term this state an "edge triplet" 
since the triplet is formed between a p orbital and the L z = —2 <i-shell orbital, unlike 
the Hund's-rule state whose triplet is formed by the two p orbitals. (The d-shell orbital, 
having greater diameter, is closer to the edge of the droplet than the p-shell.) Continuing 
the magnetic field evolution would lead eventually to the v = 1 state, followed by 
fractional occupation of the lowest-Landau-level orbitals. 

Just below 0.2 T, there is a transition between two excited states (see Fig. |TJ), 

§ Our notation | • • •) denotes appropriately antisymmetrised and normalised states. Additionally, the 
one-body state |10 j), f° r example, denotes the Fock-Darwin state with orbital quantum numbers m = 1 
and n — 0, and spin projection S z = —h/2. 
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both of which are triplets; one is the Hund's-rule state, the other the edge triplet. 
Explicitly breaking the rotational (spatial) symmetry allows mixing between different 
angular-momentum states, but not different spin states. The edge triplet and Hund's 
rule triplet are both strongly mixed near the transition. 



3. Anisotropy-induced mixing 

The nature of the anisotropy in the confinement potential of the present dots has been 
explored in Ref. where the confinement potential has been found to fit very well the 
analytic form 

2" 



T T 

Kon = 2 m ^0 



(x 2 + y 2 )+i(x- 



(4) 



The gate geometry producing such confinement is sketched in the inset of Fig. El For 
the four-electron droplet, we find 7 = 0.5 and D = 5ao, where ao = ti 2 /(m*e 2 ) is the 
effective Bohr radius. The first term in parentheses is the usual parabolic confinement, 
and the second term in parentheses is the explicit symmetry-breaking term which induces 
mixing in the angular momentum states. Because the anisotropy is of quartic order, 
non-zero matrix elements exist between states with angular momentum L z and L z ± SR, 
where 5R = 0, 1, 2, 3, 4. The respective Hilbert spaces of different spin angular momenta 
remain segregated. 

It is useful to consider the Hamiltonian as two terms 

H = H + H 1 . (5) 

H Q contains the full Coulomb interactions, the kinetic energy, and the isotropic 
(parabolic) confinement potential, whereas Hi contains only the anisotropic piece of 
the confining potential — the term proportional to 7 in (J3J. We can label the eigenstates 
of H as \i, L z , S), which denotes the i-th state with z-component of angular momentum 
L z , and spin S: 

H \i,L z ,S) = E<{ 0) (L z ,S)\i,L z ,S). (6) 

(In practise, we also use S z as a good quantum number.) Figure [T] shows the field- 
dependence of several Ef\ 

The eigenstates \j, S) of the full anisotropic Hamiltonian H in (JEJ do not preserve 
angular momentum. These states can be constructed from the parabolic states \i, L z , S): 

\j,S) = J2®t\hL z ,S). (7) 

i,L z 

The probability P(L Z ) that a given state \ j, S) has angular momentum L z is then given 
by 

P(£*) = J>M 2 , (8) 

i 

with the sum rule P(L Z ) = 1 conserving probability. 
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Our calculations of four interacting electrons in an anisotropic quantum dot, (j5J), 
and (jip. proceed through two successive exact-diagonalisation procedures. In the first 
step, we compute the energy levels and eigenstates of four interacting electrons in a 
parabolic potential, H in (jSJ). As a basis, we use antisymmetrised four-electron states 
of the 2D harmonic oscillator (Fock-Darwin) states \mna). We use the full long-range 
Coulomb interaction, whose exact matrix elements are given in Ref. We ensure the 
accuracy of our results by progressively increasing the Hilbert space and monitoring the 
convergence of the energy levels. This yields the Ef (L z , S) and the \i, L z , S) of (0). 

In the second step, we diagonalise the full anisotropic Hamiltonian. The basis in this 
step is formed from the correlated eigenstates, \i, L z , S), computed in the previous step. 
The notation u \i, L z , S)" indicates the zth eigenstate of H (which includes full Coulomb 
interactions) with angular momentum L z and spin S. Using this technique, the size of 
the Hilbert space and the degree of angular-momentum mixing can be independently 
controlled, and the convergence of the resulting energies and states can be readily 
investigated. Total spin S and z-component S z are good quantum numbers and are 
used to classify our states. 

3.1. Ground-state spin transition 

We first discuss the zero-field case. Figure ^ shows a triplet Hund's-rule ground state. 
This state has a filled s-shell and two singly-occupied p-orbitals. In contrast, the first- 
excited state is a singlet with a doubly-occupied s-orbital and one doubly-occupied 
p-orbital (L z = —1). At zero field, the singlet is degenerate with its time-reversed 
partner — a doubly-occupied p-orbital with L z = +1. Because the singlet state is 
degenerate and the triplet is not, we expect anisotropy to mix the singlet state more 
strongly than the triplet. As a consequence, the singlet-triplet gap should decrease as 
a function of anisotropy. (In other words, the lowest-energy triplet should rise faster in 
energy with anisotropy than the lowest-energy singlet.) This is indeed the shown 
in Fig. |21 In essence, the anisotropy-induced p-shell hybridisation is strong enough to 
produce two states delocalised across the L z = ±1 orbitals. The splitting between these 
hybridised states is sufficient to overcome the exchange-energy savings associated with 
the triplet state. Hund's rules are thus mitigated and the spin of the ground state can 
be tuned through the anisotropy. 

3.2. Field- dependence and angular-momentum content 

In addition to the ground-state spin, the entire field-dependence of the many-body states 
can be tuned. The states analogous to those in Fig. are shown in Fig. El The relatively 
weak field dependence of these states, is an indication of the angular-momentum mixing 
induced by the anisotropy. (Compare with Fig. which shows states with definite 
angular momentum for the same energy interval (0.8 meV).) 
The average angular momentum 

(L z ) = {j,S\L z \j,S) (9) 
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Figure 2. Main plot: Energy difference SE between the lowest-energy singlet and 
triplet states. Positive (negative) SE indicates a triplet (singlet) ground state. The 
Hund's rule ground state exists only for SE > 0. Inset: Sketch of the gate geometry 
producing anisotropic confinement, 
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Figure 3. Main plot: Filled circles indicate the two lowest-energy singlet states for 
the four-electron interacting and anisotropic droplet. Open circles are the same, but 
for triplet states. The figure is for anisotropy parameters 7 = 0.5 and D = 5ao, and 
for ujq = 1 mcV. All other parameters are taken for GaAs. Inset: Average angular 
momentum (L z ) for the lowest-energy singlet and triplet states for the same set of 
parameters, as a function of magnetic field B [T]. 



of the states shown in Fig. El are not constants but are themselves dependent on the 
magnetic field. This field dependence is shown in the inset of Fig. El At zero magnetic 
field, both states have (L z ) = 0. This is a reflection of time-reversal symmetry, 
whose consequence is that the energy of each angular-momentum state depends on 
the magnitude, not the sign, of the angular momentum. Hence, +L Z and — L z states 
are mixed equally. This degeneracy is lifted when time-reversal symmetry is explicitly 
broken, in favour of the negative angular-momentum states. 

We focus first on the singlet. In the isotropic case, the lowest-energy singlet has 
L z = —2. (See Fig. [TJ) With the anisotropy present, the inset of Fig. Elshows that (L z ) 
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Figure 4. Top: Dominant angular-momentum composition of the lowest-energy 
singlet state for the four-electron droplet. The figure is for anisotropy parameters 
7 = 0.5 and D = 5ao, and for u>o = 1 meV. the mixing occurs primarily at zero field. 
Many more angular-momentum states contribute to this state, but these all have less 
than 10% contribution. Bottom: Same plot, but for the lowest-energy triplet. In this 
case, the mixing occurs primarily at moderate (0.25 T) fields. 



drops quickly from its zero-field value of (L z ) = 0, and then, above 0.1 T, the field- 
dependence is much weaker, ranging from about (L z ) = —2 to —2.5. This indicates that 
the singlet state is primarily composed of the isotropic (L z = —2) singlet state above 
a field of approximately 0.1 T. The angular-momentum mixing is more severe below 
this field value. These results are born out by further analysis. In Fig. IH we plot the 
individual angular-momenta contribution, (J5]L for the singlet state. Only those angular 
momenta with P{L Z ) > 0.1 are shown ||. Figure 0] shows that, above 0.1 T, the isotropic 
L z = —2 state is by far the dominant contribution; the contribution of its time-reversed 
partner, L z = +2, dies out rapidly. Thus, the anisotropy does not severely mix angular 
momenta for this state above 0.1 T. 

Figure HU shows that the situation is rather different for the triplet state. Here, 
(L z ) decreases with magnetic field at a more uniform rate than does the singlet. In this 
case, significant mixing occurs between the edge triplet and the Hund's-rule state, also a 
triplet. (See Fig. [IJ) In Fig.|H we show the dominant angular momentum contributions 
to this lowest-energy triplet state, P(L Z ) > 0.1. We see that the makeup of this state 
chang GS £lS db function of magnetic field, being composed primarily of the Hund's-rule 
state at low fields, and the edge triplet at fields above about 0.25 T. Another state, 
with L z = —4, emerges above the background for fields greater than 0.3 T. Thus, the 
anisotropy induces significant angular-momentum mixing for the triplet state, although 
not at zero field. 

|| The calculations were performed by mixing 31 different angular-momentum sectors, L z — —20 to 
L z = +W. 
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4. Summary 

Well-isolated quantum dots containing a known and small number of interacting 
electrons are simple enough to be amenable to a theoretical treatment through directly 
diagonalising the relevant many-body Hamiltonian. By employing the successive exact- 
diagonalisation procedures outlined above, calculations can be extended with reasonable 
efficiency to anisotropic systems with lower symmetry. At the same time, these systems 
are amenable to experimental analysis as well. With high source-drain spectroscopy, 
excited states as well as ground states can be directly imaged in experiment. 

In the present work, we analysed ground and excited states in a four-electron 
lateral quantum dot. We find that the (many-body) long-range Coulomb repulsion 
and the (single-body) anisotropic confinement potential both play important roles and 
both must be included in any analysis. By analysing the magnetic- field dependence of 
these many-body levels, we have determined the average angular momentum of each 
level, which angular-momentum states are mixed together, how much they are mixed 
together, for both the ground and excited states of the system. 

The interplay between experiment and theory can be exploited to design gate 
geometries which in turn produce confinement potentials whose magnetic field 
dependence can be predetermined to be as desired — a flat field-dependence, for example. 
These are the elements necessary to achieve quantum-state engineering, and a necessary 
prerequisite for controlling, in real time, the coherent time evolution of individual states 
that are robust toward decohering influences. 
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